Investigation of the Substituent Effects on z-Type Pnicogen 


Bond Interaction 


XU Hui-Ying( Ew) ;CAO Sheng-Wei( #7 kb) 2; WANG Wei(E. 48) °;ZHU 


Jian-Qing( A 22) ;ZOU Jian-Wei(4h 4 2) ©; XU Xiao-Lu(V EER) PLU Yin (iti 
Jaly? 


“ (College of Biology & Environment Engineering, Zhejiang Shuren University, 
Hangzhou 310015, China;®? Zhejiang Linix Motor CO., LTD, Dongyang 322118, 


China;® Zhejiang Surveying Institute of Estuary and Coast, Hangzhou 310008, 
China;® Department of Basic Courses, Zhejiang Shuren University, Hangzhou 310015, 
China;® Ningbo Institute of Technology, Zhejiang University, Ningbo 315104, China 


ABSTRACT Intermolecular interactions between PH2Cl and Ar—R (R = H, OH, NH2, 
CH3, Br, Cl, F, CN, NO2) were calculated by using MP2/aug-cc-pVDZ quantum 
chemical method. It has been shown from our calculations that the aromatic rings with 
electron-withdrawing groups represent much weaker binding affinities than those with 
electron-donating groups. The charge-transfer interaction between PH2Cl and Ar-R 
plays an important role in the formation of pnicogen bond complexes, as revealed by 
NBO analysis. Nevertheless, AIM analysis shows that the nature of the interactions 
between PH2C1 and Ar-R is electrostatic, and the interaction energies of the complexes 
are correlated positively with the electron densities in the bond critical points (BCPs). 
RDG/ELF graphical analyses were performed to visualize the positions and strengths of 
the pnicogen bonding, as well as the spatial change of the electron localization upon the 
formation of complexes. The z-type halogen bond was also calculated, and it has been 
revealed that the z-type pnicogen bond systems are more stable than the halogen bond 
ones. 

Keywords: z-type pnicogen bonding; z-type halogen bonding; NBO; AIM; 
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1 INTRODUCTION 


Molecular interaction, closely correlated to various physicochemical properties, life 
phenomena and material structures, has been of great concern in physics, chemistry, biology, 
material science and some other fields''>!. The hydrogen bond interaction is of typical molecular 


interaction'®!, Then, the halogen bond"!, lithium bond'*! and other weak interaction were found. In 


2009, Hey-Hawkins et all proved the existence of P---P non-bond interaction via °C {'H, 3P} 
NMR experiment. Scheiner!!®'>! theoretically calculated the geometrical structure and interaction 
energy of a series of pnicogen bond systems including P---P, N---P, etc. Since then, the research on 
pnicogen bond interaction has attracted the attention of theoretical and experimental chemists!!6-20, 
Recently, Zukerman-Schepector et al!!! emphatically discussed the As---2 interaction of the 
supermolecular system, searched out 20 structures involved with the As---z interaction from the 
CSD (Cambridge Structural Database), and concluded the importance of As---z interaction in 
supermolecular construction from these characteristic structures. Based on the achievement of 
Zukerman-Schpector, et al, the interaction was calculated between ECl; (E = As, Sb, Bi) and z 
electron donors (e.g. benzene, hexafluorobenzene) by Frontera subject team!??!, and the 
significance of z-type pnicogen bond in the life science field was explored. PH3 is the simplest 
pnicogen-bond donor molecule with a weak pnicogen bond interaction. However, when one of the 
H atoms is replaced by Cl atom (PH2Cl), the P---N pnicogen bonding interaction formed between 
PH2Cl and the representative electron donor NH3 even exceeds the hydrogen bond interaction 
between water molecules!?3!, so PH2Cl is often used as the model molecule in the research on 
pnicogen bond interaction! ?51, In this paper, the interactions between PH2CI and the substituted 
benzene (Ar-R, R = H, OH, NH2, CH3, Br, Cl, F, CN, NO2) are calculated using the quantum 
chemistry method, aiming to explore the geometrical structure, electronic structure and interaction 
energy of z-type pnicogen bond between PH2C1 and the aromatic compound. Furthermore, the 
effects of the substituting group in the aromatic ring to the stability of x-type pnicogen bond 
complex are considered, in order to provide a theoretical basis for recognizing the nature of x-type 


pnicogen bonds. 


2 COMPUTATIONAL METHODS 


The structure optimizations and energy calculations of all molecules are performed by Gaussian 
03"! program package. Considering the importance of electron correlation for accurate expression 
of interaction energy, the structure optimizations of all monomers and complexes are performed by 
MP2 method with aug-cc-pVDZ basis set, and no symmetry limit is applied to the optimization. It 
has also been proved that this method is applicable in the research on pnicogen bond interaction”: 


281. The interaction energy (Ain), obtained from the complex energy deducting two monomer 


energies, is corrected to interaction energy (AE) through the BSSE (Basis set superposition 


error) correction by CP (counterpoise) procedure’?!. To further investigate the geometric and 
electronic properties of the complex systems, the second-order perturbation stabilization energy 
(AE?) and the charge transfer quantity are calculated by the NBO program in the Gaussian 03 
program package. Meanwhile, the topological properties at BCP (Bond Critical Point) are 
calculated by the program AIM 2000, and the RDG/ELF (Reduced Density Gradient and Electron 
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Localization Function) is obtained by the Multiwfn software developed by Lu, et al®°!, and used to 


analyze these systems. 


3 RESULTS AND DISCUSSION 


3.1 Geometrical structure and interaction energy 


The stable structures of the complexes combining PH2C1 with substituted aromatic compounds 
are obtained at the level MP2/aug-cc-pVDZ, and the geometrical structures are shown as Fig. 1. 
As can be seen from these geometrical structures in Fig. 1, the pnicogen bond interaction is 
generated between the P atom in PH2Cl and the z-electron in the aromatic ring. According to our 
aforesaid definitions®"!, the a angle is the included angle between the vector from the P atom to 
the hetero ring centroid and the vector at the P—Cl bond direction (see Fig. 2), and the values of 
angle a for complexes 1~9 are shown in Table 1. The angle a values of 8 complexes are almost 
180° except that of complex 3, i.e. the three points of Cl, P, and the centroid are almost at the same 
straight line. The optimization result for complex 3 shows that one of the H atoms (H17) in PH2Cl 
is also close to the aromatic ring center, i.e. a certain x-type hydrogen bond interaction is formed 
between the atom H in PH2Cl and the z-electron in the aromatic ring, so the decrease of angle a is 


caused by the interference of interaction P(15)-H(17)---z. 
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Fig. 1. Optimized geometries of complexes at the MP2/ aug-cc-pVDZ level 
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Fig. 2. Geometric model of the pnicogen bonded angle a 


The uncorrected interaction energies (AFin) and the corrected interaction energies (AES) of 


the complexes obtained at the MP2/ aug-cc-pVDZ level are listed in Table 1. Comparing the AE in 


with AEC , the differences between them are 2.7 to 3.5 kcal-mol!, so BSSE correction is 


int > > 

necessary for the x-type pnicogen bond systems. Complex 1 is the one formed by PH2C1 and the 
non-substituted benzene, whose corrected interaction energy is —3.81 kcal-mol!. The corrected 
interaction energies of complexes 2 and 3 are —7.56 and —9.65 kcal-mol', respectively, which are 
much larger than that of complex 1. As the substituted groups -OH and -NH) are electron donor 
ones, the electron density of z-electron system in the benzene ring is increased, which then leads 
to a stronger interaction between PH2Cl and z-electron, greater interaction energy of the complex 
and a more stable complex. Similarly, the electron donating capacity of the group -CH3 is weaker 
than the above two, so the interaction energy of complex 4 is smaller than that of 2 and 3. The 
corrected interaction energies of complexes 7~9 are —3.55, —3.44, and -3.45 kcal-mol"|, 
respectively, which are smaller than the interaction energy of complex 1, i.e. complexes 7~~9 are 
less stable than 1. This is because the electron-withdrawing groups, -F, -CN and -NO», reduce the 
z-electron density in the benzene ring, and result in the weakening of the interaction between P 
and z-electron. It is worth noting that the substituting groups -Br and -Cl in complexes 5 and 6 are 
electron-withdrawing groups, but their interaction energies are a little larger than that of complex 
1. The reason may be the weak hydrogen bond interaction between the hydrogen atom in PH2Cl 


and halogen atom of the substituting group. It can be clearly seen in the geometrical structure of 


Fig. 1, and also proven by the following NBO analysis. 


The halogen bond is the interaction between the electron donor and the halogen atom in the 
halogenous molecule. It was found earlier than the pnicogen bond, and the research on halogen 
bonds is more mature than the pnicogen bond. In our current research system, when the z-type 
pnicogen bond interaction can be generated between the PH2Cl molecule and the aromatic 
compound, a z-type halogen bond complex can be generated between z electron of the aromatic 
ring and atom Cl in the PH2C1 molecule. As a comparison, the molecular structures of nine typical 


halogen bond complexes are optimized at the MP2/aug-cc-pVDZ level, and then their corrected 


interaction energies (AEȘ ), corrected through the BSSE, are shown in Table 1. Comparing the 
corrected interaction energy AES of the pnicogen bond system and the corrected interaction 


energy (AES ) of the corresponding halogen bond system, it has been found that the z-type 


pnicogen bond system is more stable than the z-type halogen bond system. It is caused by the 
unbalanced distribution of positive electrostatic potential in the PH2Cl. The structures of two 
representative halogen bond complexes are shown in Fig. 3, where the positive electrostatic 
potential at the top of atom P is larger than that at the top of atom Cl in the same aromatic m 


electron donor. 


Table 1. Geometric Parameters and Interaction 


Energies of the Complexes at the MP2/aug-cc-pVDZ level?’ 


on, A an M A 
Moen -6.54 -3.81 169.1 PH,CI---Ar-H -1.58 

ee on -10.66 -7.56 169.3 PH)Cl--- Ar-OH —4.99 
E eee X = -13.09 -9.65 145.6 PH,Cl---Ar-NH> 6.00 
m nas i i -7.19 —4.28 67.1 PH,Cl--- Ar-CH3 -2.09 
pera T an C Eee ee 
ep ai ~7.22 -4.05 170.1 PH.Cl---Ar-Cl -1.97 
ao a 6.49 -3.55 71.8 PH)Cl---Ar-F -1.82 

N —6.43 -3.44 68.8 PH,Cl--- Ar-CN -2.35 
Complex? 6.63 -3.45 170.1 PH,Cl---Ar-NO> -2.74 


(H2CIP. - - Ar-NO2) 


a Interaction energies are given in kcal-mol™!. 
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Fig. 3. Molecular graphs of PH2C1::: Ar-H and PH2Cl---Ar—OH and electrostatic potential surface of 


monomers PH2Cl 
3.2 NBO analysis 


Table 2 shows the donor-acceptor orbital, the second-order perturbation stabilization energy 
(AE’), and the charge transfer quantum (Qcr) of 9 x-type pnicogen bond complexes at the level 
MP2/aug-cc-pVTZ. The second-order perturbation stabilization energy (AE?) can be obtained 


from the following equation: 


FG, iy 
E-E; 


AE’ =AEjïj=4q; 


where q; is the donor orbital occupancy, €, and €, are the diagonal elements (orbital 


energies), and F is the NBO Fock matrix element. It has been shown from the donor-acceptor 
orbital for complexes 1~9 in Table 2 that the pnicogen bonding interaction is mainly the 
interaction between the C-C z bonding orbital and the ø anti-bonding orbital of P—Cl in the 
monomer PH2C1. The three complexes 3, 5 and 6 are further involved with the hydrogen bond 
interaction in addition to the pnicogen bond interaction. The second-order perturbation 
stabilization energy values indicate that the pnicogen bond interaction, followed by the hydrogen 
bond interaction, plays a key role in the complex stability. The existence of the supportive 
hydrogen bond interaction is the fact that complexes 5 and 6 have higher interaction energies than 


complex 1 in the above analysis, and why the angle a is smaller in complex 3. 


Table 2. Natural Bond Orbital Analysis of Complexes 1~9 at the MP2/aug-cc-pVDZ Level 


(FRAT 


Complexes Donor NBOs Acceptor NBOs AE?/kcal.mol'! QOcr/e 
Complex 1 BD (2) C(1)-C(6) 3.05 
(H2CIP- -- Ar-H) BD (2) C(2)-C(3) BD* (1) P(3)-C1(16) 0.56 D pria 
Complex 2 BD (2) C(1}-C(6) 3.50 
(H2CIP--- Ar-OH) BD (2) C(2)-C(4) BD* (1) P(14)-Cl1(17) 0.62 0.01536 
Complex 3 BD (2) C(2}-C)  BD* (1) P(15)-C1(18) 4.47 anise 
(HbCIP--- Ar-NHb) BD (2) C(1)-C(6) BD* (1) H(17)-P(15) 1.02 i 


Complex 4 BD (2) C(1)-C(6) 7 2.30 

(H:CIP--- Ar-CH;) BD (2) ce) — BP* O Pt6)-Cla9) 1.10 001007, 
Complex 5 BD COC)  BD* (1) P(13}-CI(16) 2.85 Sü 

(H2CIP--- Ar-Br) Lp (3)Br(12) BD* (1) H(14)-P(13) 0.81 
Complex 6 BD COC)  BD* (1) P(13}-C1(16) 291 oain 

(H2CIP--- Ar-Cl) Lp 3)C1012) BD* (1) H(14)-P(13) 0.49 
Complex 7 BD (2) C(2)-C(6) : 3.06 

(H;CIP--- Ar-F) BD (caelo  BP* O PU3)-C1(16) 0.43 9.03077 

Complex 8 me 

Hees T BD (2) C(6)-C(7) BD* (1) P(1)-Cl(4) 2.78 0.00643 
Complex 9 BD (2) C(5)-C(6) . 2.30 

(H2CIP.... Ar-NO>) BD (2) C(7)-C(8) Eee) 0.57 9.00333 


According to the data in Table 2, the correlation between the second-order perturbation 
stabilization energy of the C-C z bonding orbital and the o anti-bonding orbital of P—Cl and the 
corrected interaction energy are mapped in Fig. 4(a), and the following relational expression is 


obtained after fitting: 
AEC = —37.0697 + 21.569 1AE? — 3.46391 (AE?) 
R=0.956, SD = 0.527, N=9 


The curve in Fig. 4(a) shows there is a bi-variable function relation between the second-order 


perturbation stabilization energy (AE?) of the C-C z bonding orbital and the c anti-bonding orbital 


CP 


it ) Other than the traditional linear relation, and 


of P-Cl and the corrected interaction energy (AE. 


it also shows the contribution of pnicogen bond interaction to the complex stability. In 
consideration of the factors of hydrogen bond interaction in the other three complexes, the 
second-order perturbation stabilization energy of hydrogen bond interaction is combined with the 


one of z-type pnicogen bond interaction, and the relationship between the combined second-order 


perturbation stabilization energy (AE’) and their corrected interaction energy (AE“") is further 


int 
mapped in Fig. 4(b). The correlation chart in Fig. 4(b) shows that the combined AEF? is linear to 
their interaction energy with the correlation coefficient (R) of 0.924, i.e. the weak interactive 
second-order perturbation stabilization energy (AE?) is positively correlated to the complex 


stability. 
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Fig. 4. Relationship between the second-order perturbation stabilization energy and interaction energy 


The molecular interaction is always accompanied by charge transfer, whose quantities in 
complexes 1~~9 are shown in Table 2. As seen in the table, the charge quantities transferred from 
the aromatic compound to PH2C1 molecule are from 6 to 16me. The charge quantity transferred is 
not so much, and it complies with the soft acid-soft alkali model (PH2C1 molecule as the soft acid 
and the aromatic compound as the soft alkali). The plot of the charge transfer quantity and the 


corrected interaction energy is shown in Fig. 5. The curvilinear equation is as below: 


AES? =-11.215+ 1852.819 Ocr -106824.298 (Qcr)? 
R=0.954, SD = 0.540, N=9 
Fig. 5 shows that the charge transfer quantities of the 9 complexes decrease with their 


interaction energies decreasing, and this result indicates that the charge transfer plays an important 


role in the complex stability. 
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Fig. 5. Relationship between charge transfer and interaction energy 


3.3 AIM analysis 


To further analyze the nature of z-type pnicogen bond interaction, the AIM (Atom in Molecule) 
theory developed by Bader is used, as the theory is often used in researches on molecular weak 
interactions'°?*>!, The typical pnicogen bond complex molecular diagram has been shown in Fig. 6, 
in which a critical point between the atom P and the aromatic ring is observed, and therefore the 
existence of x-type pnicogen bonds is proved. Table 3 shows the electron densities (pb), Laplacian 
of electron densities (Vp,), three eigenvalues (41, 42, 23) of Hessian matrix, kinetic energy 
densities (Gy), potential energy densities (Vp) and electronic energy density (Hp) of the 9 x-type 
pnicogen bond complexes at the bond critical point of pnicogen bond (BCP) at the 
MP2/aug-cc-pV DZ level. 


According to the AIM theory, charges are dispersed and the bond ionicity is stronger at BCP 
when |A; + 42| < 23 and V?p > 0, but charges are centralized and the bond covalence is stronger 
when |A; + A2| > 43 and Vp» < 0. The data in Table 3 suggest that, the Laplacian quanta (V7pp) at 
BCP are larger than zero and |A; + 22| < 43, which means the ionicity of z-type pnicogen bond is 
stronger in complexes 1~~9. The electronic energy density H» (sum of the kinetic energy density 
Gp» and the potential energy density V») is often deemed as a correct index for understanding the 
weak interaction®® 37], The interaction is a static interaction when Hy > 0, and is a covalence 
interaction when Hy < 0. All H values in Table 3 are larger than zero, which means that the 
interactions of the 9 a-type pnicogen bond complexes belong to the static interaction, and it 
complies with the above conclusion that “charges are dispersed and the bond ionicity is stronger at 
BCP when |A; + 22| < 23 and Vp» > 0”. The electron density (pp) at BCP is correlated to the bond 
strength, as the bond strength is larger if the electron density is higher. The relationship of the 
electron density (pp) at BCP and the corrected interaction energy of the complex is mapped (Fig. 


7), and the correlation coefficient R is 0.915. 


Table 3. Topological Parameters of Complexes 1~9 at the BCP Point? 


Complexes Pb Vpp A da A3 Gb Vp H, 
er 0.0103 0.0277 -0.0055  -0.0014 0.0346 0.0062 -0.0054 0.0008 
fan on 0.0121 0.0321 -0.0061  —0.0025 0.0408 0.0072 -0.0064 0.0008 
Gy arene 0.0122 0.0307 -0.0064  -0.0035 0.0406 0.0069 -0.0061 0.0008 
E 0.0105 0.0272 -0.0055  -0.0017 0.0344 0.0060 -0.0051 0.0008 
Pay am 0.0108 0.0303 -0.0054 -0.0011 0.0368 0.0067 -0.0058 0.0009 
Complex 6 0.0110 0.0305 -0.0055 -0.0014 0.0374 0.0067 -0.0058 0.0009 


(H2CIP--- Ar-Cl) 


Complex 7 


(H,CIP..Ar-F) 0.0107 0.0293 -0.0055  -0.0013 0.0361 0.0065 ~0.0057 0.0008 
Complex 8 
(HCIP.--Ar-CN) 0.0101 0.0279 -0.0051  -0.0012 0.0342 0.0061 -0.0052 0.0009 
„CIP--- 
Complex 3 0.0108 0.0298 -0.0054 -0.0017 0.0369 0.0065 -0.0055 0.0010 


(H2CIP--- Ar-NO») 


> Topological parameters are given in a.u. 
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Fig. 6. Molecular graphs of complexes 1 and 2 
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Fig. 7. Relationship between the electron density at BCPs and interaction energy 


3.4 RDG/ELF analysis 


Yang Weitao’s subject team'**! has developed a visualized method for the weak interaction 
research, through which the calculated values of reduced density gradient (RDG) function and 
sign (A2(r))p(r) of each point in the space are visualized in the RDG isosurface map. The gradient 
isosurfaces are colored according to the corresponding values of sign (A2(r))p(r), which is found to 
be a good indicator of interaction strength. We used this method in analyzing the position and 
strength of the interaction between hydrogen bonds, halogen bonds and pnicogen bonds and the 


coordination bond, and the result was promising!??*!!, The electron localization function (ELF) is 
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an important tool for electron structure researches, and it is often used to study chemical 
problems!*”] such as the molecular interaction. These two analysis methods are applied to z-type 


pnicogen bond interaction in order to visualize the interaction change. 


Fig. 8 shows the isosurface map for complexes 1 and 2 obtained from the combination of RDG 
and ELF (Fig. 8a) and the electron localization isosurface map of PH2Cl (Fig. 8b). In the case, the 
spatial position of molecular interaction can be expected, and the interaction strength is also 
observed according to the color of RDG map, in which blue represents the strong interaction, 
green the weak interaction, and red the repulsion. Comparing the color-filled RDG isosurface 
maps for complex 1 and complex 2, we can find the blue zone of complex 2 is deeper than that of 
complex 1, which suggests complex 2 has stronger pnicogen bond interaction than 1. The white 
column and vacuum ring zones in Fig. 8a represent the localization spaces of the lone-pair 
electron, z electron and valence electron. The localization space (V(P)) of the lone-pair electron on 
the atom P is obviously weakened (Fig. 8b). The reason for such phenomenon is that the 
intermolecular distance is shortened by the weak interaction, and some repulsion is generated 
between z electron and the lone-pair electron on the atom P. Such weak interaction force 
shortening the intermolecular distance is mainly the static interaction. According to the AIM 
analysis, the pnicogen bond interaction is mainly the static acting force. In other words, what is 
weakening the lone-pair electron localization space on P atom is the pnicogen bond interaction to 


a certain extent. 


Q Pe j (O) 


Complex 1 Complex 2 


Fig. 8. (a) RDG/ELF isosurface map of complexes 1 and 2; (b) ELF isosurface map (y = 0.9) of PH2C1 


4 CONCLUSION 


The geometrical structural optimization, energy calculation and topological and graphic 
analyses for the various pnicogen bond system of PH2Cl and Ar-R (R = H, OH, NH2, CH, Br, Cl, 
F, CN, NO2) have been made at level MP2/aug-cc-pVDZ. The results showed that the complex 
pnicogen bond interaction is strengthened when the substituting group in the benzene ring is the 
electron donating group, and the complex pnicogen bond interaction is weakened when the 
substituting group is the electron withdrawing group. To compare the interaction energy of the 
m-type pnicogen bond system with the halogen bond system, the interaction energy of z-type 
halogen bond system is also calculated, and it is showed that the z-type pnicogen bond system is 
more stable. NBO theory is used to analyze the correlation between the second-order perturbation 
stabilization energy and the charge transfer quantum with the interaction energy, and the result 
shows that the charge transfer plays an important role in the stability of the pnicogen bond 
complex. AIM analysis has indicated that, the nature of pnicogen bond interaction is the 
electrostatic interaction, and the electron density at BCP is positively correlated to the interaction 
energy of the pnicogen bond complex. RDG analysis showed the position and strength of 
pnicogen bond interaction, and ELF analysis indicated the change of lone-pair electron 


localization space on P atom after the complex formation. 
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RDG/ELF isosurface map can visualize the positions and strengths of the pnicogen bonding, 


as well as the spatial change of the electron localization upon the formation of complexes. 


